Pulsed beams as field probes for precision measurement 
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We describe a technique for mapping the spatial variation of static electric, static magnetic, and 
rf magnetic fields using a pulsed atomic or molecular beam. The method is demonstrated using 
a beam designed to measure the electric dipole moment of the electron. We present maps of the 
interaction region, showing sensitivity to (i) electric field variation of 1.5 V/cm at 3.3kV/cm with 
a spatial resolution of 15 mm; (ii) magnetic field variation of 5 nT with 25 mm resolution; (iii) 
radio-frequency magnetic field amplitude with 15 mm resolution. This new diagnostic technique 
is very powerful in the context of high-precision atomic and molecular physics experiments, where 
pulsed beams have not hitherto found widespread application. 

PACS numbers: 32.60.-K, 39.10.+j 
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I. INTRODUCTION 

The structure and interactions of atoms and molecules 
have been studied with great precision using atomic and 
molecular beams 111 |2j. Some examples include accu- 
rate time-keeping [3j, |j| , precise tests of physical theories 
0, Q, measurements of the constancy of physical con- 
stants 0-(l|, and searches for physics beyond the Stan- 
dard Model [Iol - [l3l ]. These kinds of experiments require 
careful control and monitoring of stray and applied fields, 
both electric and magnetic, throughout the interaction 
region of the apparatus. In this paper we develop a novel 
method for this monitoring, based on the use of pulsed 
supersonic molecular beams. Such beams are well known 
as useful spectroscopic sources for many applications in 
physics and chemistry, mainly because they can provide a 
wide variety of molecular species at low temperature and 
high average intensity [l|. Here, we point out that the 
small spatial/temporal extent of the pulses opens a way 
to map electric and magnetic fields with high spatial res- 
olution. We demonstrate this idea using a high-voltage 
parallel-plate structure that is also a radio-frequency (rf) 
transmission-line. This allows us to drive molecular hy- 
pcrfinc transitions anywhere in the interaction region. 
The line centre, lineshape and transition dynamics then 
provide information about the electric, magnetic and rf 
fields at the place where the transition occurs. 



II. THE BEAM MACHINE 

The apparatus used to demonstrate this method is 
shown schematically in Fig. [T] Its primary purpose is 
to measure the electric dipole moment of the electron 
, a very sensitive experiment that requires careful field 
mapping. A pulse of YbF molecules is produced in the 
source every 40 ms by laser-ablating a solid Yb target into 
a supersonically expanding gas jet of Ar and SFg [Til ]. 
The gas pulse passes through a 2 mm diameter skimmer, 
placed 60 mm from the source, in order to form a well- 
collimated beam. Each pulse of YbF molecules forms in a 



small region over a short time (see section IIIII) and cools 
to give a few times 10 5 molecules at the detector in the 
ro-vibronic ground state. The molecular packets have a 
mean velocity of about 590 m/s and a velocity spread of 
40 m/s FWHM, corresponding to a translational temper- 
ature of a few Kelvin. The length of a packet anywhere 
in the interaction region is therefore approximately 7% 
of its distance from the source. 



Pump 



Probe 



\ 



rfin B(»E S 



rf out 



Source skimmer Interaction region LIF 

detector 



FIG. 1: Schematic diagram of the beam machine. The molec- 
ular packet issues from the source (in the y-direction) and 
is skimmed before being optically pumped into a single hy- 
perfine state. The molecules enter the magnetically shielded 
interaction region (dashed line) and fly through a high- voltage 
capacitor that doubles as an rf transmission-line. The electric 
field points along the z-direction. The rf-magnetic field is po- 
larised along the x-direction. A small B field can be applied 
in the z-direction. After leaving the interaction region the 
molecules are detected by laser induced fluorescence. 



A single-mode continuous-wave dye laser beam, the 
pump shown in Fig.[TJ excites the molecules downstream 
from the skimmer on the A 2 Hi/2 — A 2 S + Q(0) transition 
at 552 nm. The Doppler width is suppressed by point- 
ing the pump laser beam perpendicular to the molecular 
beam. This gives a narrow enough linewidth to resolve 
the two X-state hyperfme levels F = and F = 1, sep- 
arated by 170 MHz. We tune the pump laser to excite 
the F = 1 population, which becomes depleted as a re- 
sult. Any remaining F — 1 molecules are detected at the 
end of the 1.3 m-long beam machine using a second laser 
beam of the same frequency, the probe in Fig. [TJ The 
laser-induced fluorescence (LIF) is monitored by a pho- 
tomultiplier, a typical pulse (without pump) being shown 
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in Fig. [2] versus time of arrival. The 140 ^s width of this 
pulse is essentially all due to the velocity spread of the 
beam. The temporal resolution of the LIF detection is 
7 /is, and is due principally to the 4 mm spatial width of 
the probe laser beam. 
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FIG. 2: A packet of molecules passing through the detector. 
Each molecule emits on average 1.9 photons before being lost 
to a dark state. Approximately 2% of the emitted photons 
are registered. 



The interaction region, shown in Fig.[Tl starts 450 mm 
from the skimmer and is 790 mm long. It is magnetically 
shielded to reduce the ambient field, and current-carrying 
wires inside the shields allow us to apply a magnetic field. 
Within this region there is a pair of electric field plates, 
750 mm long and 70 mm wide, with a 12 mm spacing 
accurate to better than 200 fim over the full length. 
These plates are machined from cast aluminium, then 
gold coated to improve the uniformity of the surface po- 
tential. The whole assembly is non-magnetic. With a 
field across the gap of 15 kV/cm the leakage current is 
less than 1 nA. 

The same plate structure also serves as a 34 ft trans- 
mission line for the 170 MHz radiation that drives the 
ground-state hyperfine transition. This transports the 
radiofrequency field as a fundamental TEM wave travel- 
ling parallel or antiparallel to the beam direction. Figure 
[3] shows how the rf is coupled in and out of the transmis- 
sion line using high-voltage porcelain 100 pF capacitors 
to connect to the 50 f2 non-magnetic semi-rigid cable. 
The connections are shimmed with small (~10 pF) ce- 
ramic disk capacitors to minimise the reflections at each 
end. As a packet of molecules travels along the beamline, 
we can drive a hyperfine state changing transition at any 
desired position within the plate structure by pulsing the 
rf field on for a short period of time. 



III. TIME OF FLIGHT AND SPATIAL 
RESOLUTION 

It is a key element of our mapping technique that the 
pulses of molecules have small spatial and temporal ex- 
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FIG. 3: rf transmission line coupling. Both ends of the capac- 
itor plates are coupled to 50 fl co-axial feeds in this manner, 
rf is injected in one end of the capacitor and coupled out of the 
other. The direction can be reversed to check for systematic 
effects. 



tent at the source. A very similar source used to make 
LiH produces the molecules in the 5 fis following the firing 
of the ablation laser. This was determined from studies 
with dual ablation pulses having a variable time separa- 
tion [TB|. We believe that the YbF molecules are formed 
within a similar temporal window. The initial size of the 
YbF pulse is measured by comparing the temporal profile 
of the molecular packet at two distances from the source, 
using the pump beam as one LIF detector and the probe 
as the other. Extrapolation back to the source then gives 
the initial spatial extent as 11 mm for an initial 5 fis du- 
ration and, in any case, not more than 12 mm. 

In demonstrating our field mapping method, we have 
typically used rf pulses of 18 /is duration. During this 
short time the molecules move only 11 mm, allowing the 
rf transition to provide almost a freeze-frame picture of 
the molecular pulse at a particular position inside the in- 
teraction region. By studying the centre-frequency, line- 
shape and power dependence of the hyperfine transition 
we can infer the electric, magnetic and rf fields at that po- 
sition. Repeating the experiment for a range of positions 
allows us to build a map of the field in the interaction 
region (l8j . 

By binning the molecules according to arrival time, we 
can obtain maps with even finer spatial resolution. In 
the limit where all of the molecules are created at the 
same time and in the same place, the arrival time, t arr , 
is simply a measure of the molecule's velocity. From this 
we can infer its position z' r j at the time t r f of the rf pulse: 



-rf 



tc. 



where l m is the length of the machine from source to 
detector. The spatial resolution is limited by three fac- 
tors: the detector has an active volume of non-negligible 
size and hence limited time resolution, the rf pulses are 
of finite length, and the beam source is not point-like in 
space or time. The uncertainty from the initial extent of 
the molecular pulse can be estimated as follows. Assume 
a molecule is formed at time At and position Al with 
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velocity v. The time of arrival at the detector will be 



to 



a; 



At 



It is straightforward to show that our estimate of position 
at time t r f will differ from the molecule's true position, 

z rf , by 



Az r f — 



'rf 



t arr trf L m At t arr Al 



t, 



tnrr At 



Assuming At and Al to be uncorrelated, with standard 
deviations at and ai respectively and with zero mean, the 
uncertainty in the estimate of the position is given by 
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(1) 



For our source, Eq. [T] yields an uncertainty of 11 mm, 
which is equal to that from the rf pulse duration, giving 
a net positional uncertainty of 15 mm (the detector reso- 
lution is small). We note that the effect of finite detector 
size and rf pulse duration could be removed by an appro- 
priate deconvolution, but that approach is not pursued 
here. 



IV. FIELD MAPPING 

A. Electric field measurement 

After the F = 1 state population has been depleted 
by passing through the probe laser, it is replenished by 
the application of an rf pulse, which transfers population 
from the F = state. This causes an increase in the 
probe fluorescence, as shown in the rf frequency scan of 
Fig. [4] When considering the state of the molecule, it 
is appropriate to chose the static electric field direction 
as the axis of quantization because the electric dipole in- 
teraction is by far the strongest coupling to the external 
fields. Being a TEM wave, the rf magnetic field between 
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FIG. 4: Spectrum of the ground-state hyperfine transition. 



the plates is accurately perpendicular to the static elec- 
tric field. As a result of this polarization, the rf field 
drives the transitions (F = 0) — > (F = l,mj? = ±1). 
Transitions to (F = 1, nip = 0) are strongly suppressed 
because of the polarisation and are in any case Stark 
shifted to a different resonant frequency. For the reso- 
nance shown in Fig. [4j the rf pulse was applied when the 
molecules were near the middle of the interaction region 
and the pulse amplitude was chosen to give the molecules 
a 7r-flip on the Bloch sphere at resonance. The transition 
probability peaks at 170.821 MHz, not at the field-free 
frequency of 170.254 MHz. Because the Stark effect of 
the YbF hyperfine levels is accurately known [lj|, this 
shift allows us to determine the strength of the electric 
field at the place where the molecules are. The Doppler 
shift is only a few hundred Hertz and is therefore insignif- 
icant at this level of accuracy. The curve shown in Fig. 
[4] is the lineshape calculated for a square rf pulse. This 
agrees quite well with the data points, though not per- 
fectly. The deviations, mainly in the wings, are caused by 
the departure from a perfectly square rf pulse profile due 
to the few-microsecond settling time of the rf switches. 

In figured! the LIF signal is the total photons counted 
per beam pulse, i.e. integrated over arrival time. Figure 



5(a) shows the same data, but here the LIF signal is plot- 
It can be 



ted versus arrival time as well as frequency [H 
seen that the lineshape and line-centre are essentially in- 
dependent of arrival time, indicating that molecules at 
the front of the packet (early arrivals) experience the 
same static electric field during the rf pulse as those at 
the back (late arrivals). This shows that the electric field 
is homogeneous in the centre of the interaction region. 
In Fig. 
leaves t 



5(b) the rf pulse is applied just as the molecule 



he high-voltage plates and the frequency scan is 
zoomed in on the central fringe of the resonance. Here 
we see that the centre frequency (indicated by the dots) 
is constant for the slow molecules, which are still between 
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FIG. 5: (a) Spectrum of the ground state hyperfine transi- 
tion as a function of frequency and arrival time. Integrating 
this data over arrival time yields Fig. [4] Integrating along the 
frequency axis gives a time-of-flight profile, similar to that of 
Fig- E (b) Hyperfine transition near the end of the interac- 
tion region. Note that both axes are expanded and shifted 
compared to those in (a). Circles show the fitted transition 
frequencies for several time slices. 
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the plates when the rf is pulsed on, but shifts dramati- 
cally for the fast molecules which explore the region of 
decreasing field at the exit of the plates. 

This method of binning the data by arrival time allows 
us to map out a region, typically 6 cm in length, with a 
single frequency scan. In order to map the electric field 
over the whole interaction region, we repeat the scan for 
several rf pulse delays to make a sequence of short maps 
that can be stitched together, as shown in Fig. [5J This 
data set was collected in approximately 15 minutes. Af- 
ter fitting the line centre at each position and converting 
this to electric field [l6| , we obtain the field map plotted 
in Fig. [7] This shows very clearly a systematic ±110 /im 
variation of the plate spacing, which we can correlate 
with the design of the support mechanism. The impor- 
tance of this measurement is that it shows the variation 
of electric field directly and in situ, which mechanical 
measurements on the bench cannot do. 

The hyperfine transition frequency is sensitive to the 
electric field up to very high strengths: for YbF the 
gradient is 100 Hz per V/cm at 30 kV/cm. Because 
of this, polar molecules such as YbF offer unique sen- 
sitivity for probing very high fields. Uncertainty in the 
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FIG. 7: Electric field as a function of position. The applied 
potential was roughly ±2 kV over the 12 mm gap. 



electric dipole moment of the molecule limits the abso- 
lute accuracy to approximately 1% in the case of YbF 
16;] . However, small relative field measurements are typi- 
cally limited only by the accuracy of determining changes 
in the centre of the line. This in turn depends on the 
linewidth associated with the duration of the rf pulse. 
In the present experiment, using scans with 15 mm spa- 
tial resolution (18/Lts pulse), we achieve an uncertainty 
of 500 Hz with 30s of measurement, corresponding to a 
relative field uncertainty of only 1.5 V/cm at 3.3 kV/cm. 



B. Static magnetic field 

If there is a magnetic field component B z parallel to 
the electric field, the (F = 1, mp = ±1) energy levels are 
split by the Zeeman interaction. Although this does not 
shift the centre of the resonance, it does cause a splitting 
of the line, which can be used to measure B z . By con- 
trast, a perpendicular component of magnetic field has 
no significant effect on the lineshape and produces only 
a very small shift of the centre. This insensitivity to the 
perpendicular field is due to the strong tensor Stark ef- 
fect, which shifts the (F = 1, mp = 0) state far away in 
energy compared with the off-diagonal magnetic coupling 
fill ] . By switching off the electric field it is also possible 
to measure B y , but we do not demonstrate that here. 

In order to demonstrate magnetic field mapping, we 
apply a small static magnetic field, roughly parallel to 
the strong electric field, by passing a current / through 
through copper wires that run parallel to the beam in- 
side the innermost magnetic shield. In addition there 
is a residual ambient field from imperfect shielding, and 
perhaps from ferromagnetic impurities within the inter- 
action region. 

Figure [5] shows the resonance line measured near the 
centre of the machine with a Zeeman splitting of 25 kHz 
due to the application of a 0.9 ^T field. The rf pulse du- 
ration was lengthened to 50 fxs so that this splitting could 
be clearly resolved. The theoretical line profile describes 
the central part of the line quite adequately, but does not 
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FIG. 8: Zeeman split hyperfine transition. The fit is to a sum 
of two square-pulse transition lineshapes. 



fit well in the wings because it is once again based on a 
square rf pulse. In order to separate the ambient field 
from the applied field, we vary the current I and fit the 
field to the form B z = (31 + Bo. Once this is completed 
at a given position, the rf pulse timing can be changed 
and the measurements repeated to build up maps of /3 
and Bq along the beamline. 

Figure [§] shows the results obtained after approxi- 
mately 4 hours of data collection. The Zeeman splitting 
is fitted typically to within 100 Hz. This fitting error 
leads to an uncertainty in Bq of about 5 nT, and an un- 
certainty in (3 of 100 pT/mA. In addition, the absolute 
calibration of the field is limited by a 2% uncertainty 
in the g-factor. In Fig. 9(a) fi decreases at each end 



of the interaction region because the field rotates away 
from the z-axis to satisfy the boundary condition on the 
end-caps of the magnetic shield. The background field 
shown in Fig. |9(b)| exhibits a sharp increase in magni- 
tude at both ends of the interaction region. This was not 
expected and requires further investigation, though we 
surmise that it is due to inadequate demagnetisation of 
the magnetic shield end caps. 



C. rf magnetic field 

In the previous two sections, the centre and lineshape 
of the rf transition have informed us about the static 
electric and magnetic fields. In this section we use the 
dynamical evolution of the transition to provide informa- 
tion about the rf magnetic field itself. 

LIF detector signal (photons / pulse) 
1600 

1400 *• '.- 

1200 "' 
1000 »' 
800 

600 ,* 
t 

400 r \ 

12 3 4 

rf amplitude (arb.) 

FIG. 10: Number of molecules in the F = 1 state versus rf 
amplitude. 



The rf field supported by the parallel plate transmis- 
sion line could be a pure travelling wave in principle, 
but imperfect impedance matching at the ends inevitably 
produces some standing- wave component. It is therefore 
of interest to map the rf field strength as a function of 
position along the beamline. For this measurement, we 
fix the rf frequency at resonance and scan the amplitude 
to observe the Rabi flopping of the molecular population 
between the F = and F = 1 hyperfine states. Figure 
[TTJl shows such a scan covering almost An of Rabi flop- 
ping, with the rf pulse timed to excite molecules at the 
centre of the interaction region. The very good visibil- 
ity of these fringes indicates that the rf field is rather 
homogeneous over the beam. In Fig. 11(a) we explore 
the longitudinal direction by replotting the data versus 
arrival time as well as field amplitude, much as we did 




rf amplitude (arb.) rf amplitude (arb.) 



FIG. 9: Maps of (a) applied field coefficient j3 and (b) ambient ( a ) Homogeneous field (b) Inhomogeneous field 

field Bo in the interaction region. 

FIG. 11: Rabi-fiopping, resolved by arrival time: (a) driven 
by the homogeneous field of the transmission line, (b) driven 
by the inhomogeneous field of an 8 cm diameter loop. 
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make a single 7r-flip. The variation of the rf magnetic field 
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FIG. 12: (a) Points: magnetic field of rf loop derived from 
data of fig. 11(b) Line: fit to dc field of a loop. Peak strength 
is the only parameter, (b) Points: rf field in the transmission 
line. Line: sine-curve corresponding to the field amplitude for 
a 170.820 MHz standing wave. The amplitude and phase of 
the standing wave, and amplitude of the travelling wave were 
floated in the fit. Note the different x-axis scale compared to 
(a). 



amplitude reconstructed from this is shown in Fig. 12(a) 



The line superimposed on the data is the simple magnetic 
field variation expected for a DC field on the axis of a 
loop. 

Returning to the transmission line, we have measured 
the rf field along the whole length of the interaction re- 
gion by varying the timing of the pulse and making many 
overlapping scans. Figure |12(b)| shows the result and 
clearly reveals the standing-wave component of the field. 
Superimposed on the data is a sine curve of half the free- 
space wavelength (87.5 cm), which fits well to the ob- 
served modulation. From this map, we estimate that the 
reflected amplitude at each end of the transmission line 
is 20%, as one would expect in view of the impedance 
step between 50 fl and 34 51. 



in Fig. 
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This shows that molecules in different posi- 
ong the length of the packet have the same Rabi 
flopping rate, further demonstrating the uniformity of 
the field. 



In contrast to this, Fig. 11(b) shows the result obtained 



with an inhomogeneous rf field generated by a copper 
loop encircling the field plates. The rf pulse timing was 
chosen so as to centre the molecules in the plane of the 
loop. One sees that molecules at the centre of the packet, 
with an arrival time of approximately 2.24ms, undergo a 
Rabi rotation of 57r at the highest rf amplitude plotted, 
whereas molecules at the edges of the packet can barely 



V. CONCLUSION 

We have described a set of sophisticated diagnostic 
techniques that can be applied to pulsed atomic and 
molecular beams. These techniques have been developed 
in our lab to control and measure systematic effects as- 
sociated with our electron-edm measurement. However, 
the main point of this article is to note the wider appli- 
cability of these methods to any experiment using pulsed 
atomic or molecular beams. In particular, these methods 
are very powerful in the context of high-precision atomic 
and molecular physics experiments, where pulsed beams 
have not hitherto found widespread application. 
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